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Abstract 

fvl .We study (analytic) finite-size corrections in the dense polymer model on the strip by perturbing the critical Hamil- 
y—i 'tonian with irrelevant operators belonging to the tower of the identity. We generalize the perturbation expansion to 
^D .include Jordan cells, and examine whether the finite-size corrections are sensitive to the properties of indecomposable 
^^ 'representations appearing in the conformal spectrum, in particular their indecomposability parameters. We find, at first 



\0' 



.order, that the corrections do not depend on these parameters nor even on the presence of Jordan cells. Though the 
^corrections themselves are not universal, the ratios are universal and correctly reproduced by the conformal perturbative 



.approach, to first order. 
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.1. Introduction 

] Universality, scaling, and finite-size corrections for crit- 
ical systems have attracted much attention in recent deca- 
des [ij, [3, i, H i, @, 0, il ■ Although many theoretical results 
■are now known about the critical exponents and universal 
relations among the leading critical amplitudes, not much 
information is available on ratios among the amplitudes in 
^finite-size correction terms [3|. New universal amplitude 
'ratios have been recently presented for the Ising model 
|[6|. In this paper, we compute finite-size corrections in 
ithe critical dense polymer model [9|, and compare them 
with the results obtained in perturbation theory. Com- 
ipared to previous analyses, the new feature here is the 
^presence of indecomposable representations in the confor- 
imal spectrum [9|, |lO|, |ll| , which first requires to adapt the 
'perturbation theory to include such representations. In 
^particular it raises the question as to whether the finite- 
'size corrections depend on their indecomposability param- 
icters 



SEQEE 



Conformal invariance implies that on an infinitely long 
strip of finite width N, the eigenstates of the critical trans- 
fer matrix, associated to conformal states, have energies 
which scale with N like [17| 



EniN)^Nfb + fs + 



^C(A„ - c/24) 

N 



+ 0{N-^), (1) 



where the bulk free energy density /(,, the surface free en- 
ergy fs , and the anisotropy factor ^ are non- universal con- 
stants; in contrast, the central charge c and the weights 
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A„ of the conformal eigenstates, are universal (although 
boundary dependent for A„). 

Higher order correction terms in ([T]) are nonuniversal; 
however it has been suggested recently IB], and partially 
checked for the Ising universality class (c — 1/2), that the 
asymptotic expansion of the critical eigenvalues has the 
form (n = will refer to the groundstate) 



EniN) = Nfi, + fs + J2 



An) 



N^P- 



(2) 



and that the amplitude ratios aj," /flp are universal and 
depend only on the boundary conditions [y, |l8| . The case 
p — 1 readily follows from ([T]). 

For critical dense polymers, it is not difficult to see from 
the results of i9| that the finite-size energy levels have an 
asymptotic expansion like in ^. The general calculation 
of all coefficients is a tedious exercise, which can however 
be completed for an infinite number of levels (among these, 
the lowest-lying levels, as well as the lowest-lying logarith- 
mic pairs, in each representation). Moreover we check that 
the first coefficients Oj match the perturbative values, 
and that the ratios 0/2/112 ^'"^ universal. 



2. Critical dense polymers 

We consider the model of critical dense polymers on a 
square grid, as formulated in [9|. The degrees of freedom 
are localized on the elementary faces of the lattice, which 
can be in two different configurations, see Eq. (2.1) of [9|. 
The model on a strip, of width N and height 2M, has a 
double-row transfer matrix D(u). The partition function 
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is given by 



^N.M 



= TrD(«)^-^=^e-2A^^"(^^"), 



(3) 



where the sum is over ah eigenvalues of D(u), written 
u is a spectral parameter related to the 



as e-2S„(JV;«). 

anisotropy factor by ^ = sin 2u. 

The configuration space on which the transfer matrix 
acts can be divided into sectors Cnj, labelled by an integer 
^ > 0, of the same parity as N. The integer £ is the number 
of defect lines, and can be seen as fixing a specific boundary 
condition. The range of £ is finite when A^ is finite, but will 
be considered as unbounded in the limit A^ — >■ oo. It was 
found in [9[ that the transfer matrix is diagonalizable when 
restricted to each sector separately, although analytic and 
numerical analyses indicated that it is not in the full space 
when N is even. This was confirmed in |10l . Il9j , where it 
was shown that the sectors i and i + 2, for £ = mod 4, 
are linked up by rank 2 Jordan blocks, while there are no 
Jordan cells when N is odd. 

All eigenvalues have been determined in [9], for any 
finite value of N. They are given by 



Er, 



N-2 



E i°g (1+^^. 



j = N mod 2 



: sm ) ( 1 + C Ml sm , 

2N'^ ^^^ 2N' 



(4) 

where the summation includes the integers j with the same 
parity as N. The eigenvalues also depend on parameters 
£j, fij equal to ±1, although not all choices of Ej, /ij = ±1 
are allowed. Exactly which sequences of +1,-1 corre- 
spond to actual eigenvalues, and for which sector, are given 
by the selection rules conjectured in [9|, and proved very 
recently in [20J. What we need is summarized below. 

A main result from [9] is that the set of eigenvalues 
of the transfer matrix in the sector CN,e is such that it 
leads, through ([T]), to a set of conformal weights A„ whose 
values and degeneracies exactly match those of a highest 



_ (£-2)! 



with 



weight representation Vi,s of weight hi^s 
s = £+1. As a representation of the Virasoro algebra, Vi,s 
is believed to be the quotient of the highest weight Verma 
module Vi^s by the singular vector at level s [19[. 

The overall groundstate Eg, found in Cn.o or Cn^i, 
corresponds to all Sj = Hj = 1 in ^. From the Euler- 
MacLaurin summation formula, the asymptotic expansion 
of Eq takes the form ^ with the coefficients 



a. 



(0) 



i2p)l 



A 



■2p~l 



(C), 



N 
a — — mod 1, 



(5) 



where B„(z) are the Bernoulli polynomials, a = 0,^ and 
A2p-i(C) s-re polynomials of degree 2p — 1, obtained from 
the following expansion. 



log(l + Csinx) 



yXniC), 

n=l 



n\ 



(6) 



The first few odd order polynomials read Ai(0 = C, A3(C) = 
2C^ - C, hiO = 24C'^ - 20C^ + C, see also [Ilj. 

For p = 1 in particular, one finds Ap — c/24 = 1/12 
for N even and —1/24 for N odd. Assuming c = — 2, this 
gives Aq = /ii.i = for A^ even, and Aq = /11.2 = —1/8 
for N odd. 

For p = 2 one finds from the previous formula 



4"^ = Y(A^-Tio)A3(C). 



(7) 



The excited levels are obtained by switching the Sj , fij 
from +1 to —1 in a way allowed by the selection rules 
[9|, |20| . The specific levels we will consider here are the two 
lowest- lying levels in each sector Cn,i, denoted Eg and 
E[ , for all ^ > 1 except £ = 2 (which is somewhat special 
and requires a separate treatment; the following checks 
however work for them too). They are non-degenerate 
(within their sector) and are obtained by setting to —1 the 
parameters fij with the following indices (keeping those of 
the appropriate parity). 



l<j< 



for E, 



' 



l<j ^e-'2<£ ioT E\ 



W 



(8) 
(9) 



Working out the asymptotic expansion of Er for 



(e) 



0, 1, yields the first coefficients (with no assumption on c) 



and 



ai 



(i?W)=7r(/ii,,-(-r-f j^)Ai(C) 



a2K0 = y(/^L-Tio)A3(C), 



(10) 



(11) 



^M'') - y (^l 



6/ii,,-fii§)A3(C). (12) 



Before comparing them with the perturbative approach, it 
is useful to know where the Jordan blocks appear when one 
considers the full configuration space, namely the union of 
all sectors. 

As mentioned above, Jordan cells appear only for N 
even, in which case the full configuration space is given by 
^i>o,evenJ~'N,e- It tums out that within each sector, the 
eigenstates are generically non-degenerate, implying that 
the transfer matrix is diagonalizable in each sector sepa- 
rately; however every sector Cn,i has a subset of eigen- 
states which are pairwise degenerate with eigenstates of 
t~-N,i-2, while the complementary subset is pairwise degen- 
erate with eigenstates of Cnj+2', these double degeneracies 
occur for all values of N (even) and C. 

In the conformal limit (the order N~^ in the expansion 
of the eigenvalues), the degeneracy of the eigenstates in- 
creases. For N odd, each sector i odd carries an irreducible 
highest weight representation Vi_£+i. For N even, it is be- 
lieved [9|, |lO|, |l9[ that the pairs of sectors £ and £ + 2, for 
£ = mod 4, and their associated representations Vi^^+i 
and Vi^i+3, are tied tog eth er by rank 2 Jordan cells to form 
staggered modules 1^ [isl, [ij] . 




Figure 1: (Color online) The cone-like figures represent the represen- 
tations Vi m-i (left) and Vi ^+3 (right). The numbers on the right 
show the levels at which the various fields occur. 



Fig. 1 schematically represents the two representations 
Vi,£-i-i and V1J+3. The left one has highest weight state ^, 
a first singular vector 0, and a second singular vector 0', 
set to zero, and similarly on the right for f/;, p and p' , the 
latter being also set to zero. The states in the two blue 
regions are conformal states coming from eigenstates of 
the transfer matrix which are pairwise degenerate for all 
(large enough) values of N; they form logarithmic pairs 
and are related by the non-diagonalizable action of Lq, 
like f.i. (Lo - A)</> = 0, (Lo - A)^ == 0. 

The full projective representation, denoted by 7?.i+f /2, 
spans the two sectors £-n,i and Cn,i+2, and can be parti- 
tioned into two disjoint subsets : 5*1 contains eigenstates of 
Lq which are not parts of a logarithmic pair (green and red 
regions), and ^2 contains logarithmic pairs (in blue). Each 
such representation 7^i+f/2 comes with an intrinsic param- 



eter /3i_|.£/2j which fixes its equivalence class |12l. ll3l . ll4ll5l |. 



3. Conformal perturbation theory 

The finite-size corrections in ([2]) can in principle be 
computed in perturbation theory [22j . In general, the crit- 
ical lattice Hamiltonian will contain correction terms to 
the fixed-point Hamiltonian He = -^(-^o ^ jz)^ 



H = CH,,+ y ^ Sk / Mv)dv 



k 



N 



Perturbation techniques have been studied for a long 
time, and successfully applied in concrete models [2J| . Wri- 
te H = Ho+ gV. A basic result states that the first-order 
correction to the energy Ec^n of a non-degenerate eigen- 
state \n) of Hq is given by {n\gV\n). Here it is sufficient 
that the state \n) be non-degenerate within its own repre- 
sentation since (p preserves each representation. 

This scheme can be readily applied to the critical dense 
polymer model when N is odd, or to the two lowest-lying 
states in the indecomposable representations TZk (C and 
L_i^) when N is even, but should be adapted for states in 
Jordan cells. In addition the representations TZk carry an 
intrinsic parameter j3k, which may affect the corrections. 

It has been appropriately emphasized in Ig, [ll| that 
the conformal states should be equipped with a scalar 
product for which the conformal Hamiltonian is hermi- 
tian. When the Hamiltonian is not diagonalizable, this 
bilinear form cannot be positive definite and differs from 
the canonical scalar product usually used in the lattice 
configuration space. For a typical rank 2 Jordan cell in 
the subset 5*2 cTZk, one can take the scalar product to be 
such that {(j)\(j)) = 0, (0|V') = (^10) = /3j (^IV') = 7 where 
(3 = /3fc if {(f>,ip) is the lowest logarithmic pair of TZk, oth- 
erwise /3 is a calculable multiple of (3k for a descendant 
logarithmic pair |llj]. For states in Si, the usual scalar 
product is used. 

Assume that 0, ?/; are the conformal states of a non- 
degenerate logarithmic pair (only pair at their level). 



{Ho - Xo)4> = 0, {Ho - Ao)V = 0, 



(14) 



with Aq = A being their conformal weight (we take Hq = 
Lo for simplicity). Since (j),^p are the conformal shadows 
of states which are degenerate for all (sufficiently large) 
values of N, we may assume that the Jordan cells are pre- 
served by the perturbation, to any order, that is. 



(i7-A)$ = 0, (i/-A)* = $. 



(15) 



Expanding in power series in g, X 

(13) + g(?!)(^^ +... and * = V' + ffV''^^ - 
order the consistency condition (0| 



= A + 5A1 -f ..., $ = 
f . . ., we obtain at first 
V\(h) = and 



where the ipk are irrelevant quasi-primary fields, and gk 
are non-universal constants. On general grounds, a per- 
turbation ifk with scaling dimension Xk brings a first-order 
correction to the energy proportional to N^~^''. 

Fields in the tower of the identity are available in any 
conformal theory. On the strip, the lowest contributing 
fielqll is ip — £^2(^)1 of dimension 4 [23[; it corrects the 
conformal energies by a N^^ term, and should reproduce 



the first correction coefficient 



An) 



in (0). 



iThe 



Li^iv) - iL-iiv) 



combination of dimension 4 is 
however the terms L—„{v) for n > 3 give no 



correct quasi-primary 
3 



contribution to any order in perturbation [23 



/3Ai 



\V\i;) = imi^)- 



(16) 



Since ijj is the only state to have a non-zero scalar product 
with (j>, the r.h.s. is proportional to /?, and so Ai does not 
depend on /3. 

The required matrix element can be computed by the 
procedure used in [23| for highest weight representations. 
To avoid confusions, we denote hereafter by ((/>fe,'0fe) the 
lowest logarithmic pair in TZk, with conformal weight h. 
Thus we view the pair (0, ip) as descendants fields at level 
r of {(j)k,i'k), A = h + r. 

The main idea is to compare the chiral 3-point function 
on a strip of width N, 



F^kV^Pki^i) = {ipk{wi) L^_2{w2)^k{wz)) 



strip 



(17) 



with its expression in the operator formahsm. The pre- 
vious correlator is first computed on the upper-half plane 
and then conformally transformed to the strip. It can be 
shown that the logarithmic part of F^^^^^ is 



F, 



'■^kv^A^i) = Pk\ log— - 21og(l -] [f^ip^(zi), 

log ^ Zl V Z\ / ) 

(18) 
where F^i^^r is the strip correlator ((TT)) with ipk replaced 
by a primary field tt of the same weight h as V'/ci and 
where Zi = exp {irWi/N) relates the planar variables to 
the strip coordinates. The calculation of Ft^^t^ is standard 
but requires the transformation law of the descendant field 
L'^_2iw). Using the results of [25[ on the transformation 
laws of general descendant fields, one finds the following 
explicit function of x = Z2I ~^\ and y — z-^j zi^ 






(1 - .Ty)2'i 

jl-xy)^ 
(l-a;)4(l-j/)4- 



3h 



{l~xy)^{l + xy) 



10 



{l~xyy 



(l-x)3(l 
c 5c + 22 



12 fl 



X)2(l-y)2 

In the operator formalism, and for w 
correlator is 



24 120 



}. (19) 
+ iv, the 



F^,^^, = (0|V'fc(0,«i)T"i-"^^(0,«2)T"^-"^^fc(0,z;3)|0), 

(20) 
with the transfer matrix T = e ^^o/n ^ Inserting a com- 
plete set of states of TZk , the contribution coming from the 
pair ((/), -0) reads 

i^^(0|^fc(0,i>i)r"i-"^P^(0,7;2)Pr"=-"^^fc(0,z;3)|0), 

(21) 
where 

|V')(0| ^ |<^)(^| J^){<f>\ ,„„. 

is a projector on the subspace spanned by 4>,'>p, namely 
P\tjj) = 10), P|(/,) = 10) and P|n) = for any other state 
\n). Using the action of the transfer operator, namely 
e^io|0) = e^^|0) and e^-^°|0^) = e^'^(IV') + x\<p)), we find 
that (|21l) has a term quadratic in the w's, 
^2 



p--^^(«l-«3) 



(Wl -U2)(U2 -Us) 



/327V2 

x(O|^fc(O,«i)|0) (01^(0, 1^2)10) (0|0fe(O,z;3)|O). (23) 

As u '-^ logz, this term would require a corresponding 
log term in the 3-point function F^,,^^,,, which does not 
exist. We thus have (0|i^(O, U2)|0) = 0, which fulfills the 
consistency condition we found in perturbation theory. 
Looking now at the terms linear in the u's, we find 

^ p--grA(Mi-ti3) y 



P'^N 



{(ui-«2)(O|^fc(O,Vl)|0)(0|<^(O,«2)|l^)(0|^fe(O,W3)|O) 
+ (u2-li3)(O|V^fc(O,l>i)|0)(V'|(^(O,«2)|0)(0|^fc(O, 1^3)10)}. 

(24) 



The dependence in U2 must vanish, which implies the equal- 
ity (0|i^(O, W2)1'0) — ('0|<^(O, U2)|0), SO that a single term 
remains. 



/327V 



{o\Mo,yiM) (01(^(0, «2)|V') (0|^fc(o,«3)|o). (25) 

Since e""^^"^^"^-* ~ z^/zi, a comparison with (jlSp implies, 
provided the pair (0, -0) is not degenerate, that the prod- 
uct of the three matrix elements in the previous equation 
equals the coefficient of (z^/zi)^ = (xy)^ in (i'^tikF-n^Tr, 
and does not depend on V2. 

Repeating the same steps for the 2-point function 

F^^^Awi,W'i) = (■0fc(wi) V'fc(w3))strip, (26) 

we obtain similarly that its logarithmic part is given by 



log 



/3fc{ log £^ - 21og (l - 'f) }f^.(z,), 

(27) 

where Ftttt is the strip 2-point function for a primary field 
of weight h, 



i'-KTT\Zi) ~ ( Af J 



2h 



{zi-z^Y^^' 



(28) 



a function of the single argument xy. 

Comparing this function with its corresponding expres- 
sion in the operator formalism, we can conclude that the 
product (01-0^(0, wi) 1 0) (0|V'fc(O,W3)|O) is equal to the coef- 
ficient of (z3/zi)^ in (3(3kF^^. 

Combining the above two results on the triple and dou- 
ble products, we obtain that their ratio yields 



(0|(^(O,t;2)|V') = (01^(0,0) I V)=/3- 



^7r^7T\X, y ) 



{xy)^ 



F-jririxy) 



(29) 



{xy)^ 



This formula shows that the first correction Ai to the con- 
formal energy of ^0 can be computed as if tp was a (non- 
degenerate) descendant |A) of a primary field tt, 



Ai 






N- 



1 



N 



^TTTT (xy) 



dv (0|i^(O, v)\il>) 
= (A|V|A). 



(30) 



{xy}^ 



Being degenerate with ^, the first correction for is equal 
to that of ip- Thus all Jordan couplings between sectors 
for N even can be neglected for the calculation of Ai, and 
each sector Cn.i can be considered separately, as carrying a 
highest weight representation Vi.s- Note that the formula 
(|29)1 requires hy^O. 

At higher orders in perturbation, we have no general 
argument to show that the same conclusion holds. On 



the other hand, we expect that the first order corrections 
coming from more irrelevant perturbations in the tower of 
the identity show the same independence to the Jordan 
structure. 

So let I A) — \h, r) be a non-degenerate, level r descen- 
dant state in a representation of highest weight h. The nor- 
malizeco matrix elements Cai^a = (-^/7'')''(A|(p(0, 0)|A) 
for tf — {L^i^i) can be computed from (|29p by using the 
explicit functions -FjripTr and F-^^^ above, and reproduce, for 
h ^ 0, the results first found in [231], 
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CavA = ih + r)(^h + 



r{2h + r){5h + 1) 



{h + l){2h + l) 
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- (5c- 



(31) 



We obtain from this the first correction to the conformal 
energy of a state |A) as 



ffAi = {A\gV\A) = ^Ca^a- 



(32) 



W 



We may now apply this to the two lowest levels Eq 

and El in each sector jCN,e (the only ones to be gener- 
ically non-degenerate). For those states, h = hi^g and 
r = or r = 1, so that we find from the formula ([5T|) with 
c=-2. 



1 

120 



\ 9'^ (u-2. 



r = 1. 



(33) 
(34) 



The corrections obtained earlier in 0, pT|) and (IT^ 
from the lattice transfer matrix are clearly reproduced 
by setting g = A3(C)/67r. We note that c = — 2 is the 
only value for which the perturbation theory correctly re- 
produces the lattice finite-size corrections. If we do not 
fix c, we see from ([TUl) that the weight h is shifted to 
h = his + 24 ' ^^'^ from the formula (J31L the first or- 
der Ai for r — and r — 1 gets an extra term, inside the 
parentheses, equal to -^440 ^^-^ — j^^^^q ' respectively. 
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^The normalization in |23|1 is {N/2it)'^ because the calculation 
is performed on a cylinder, in which case the planar and cylinder 
variables are related by 2 = exp {2nw/N). 



